Accepted for publication in ApJ - February 21, 2008 



Morphological evidence for azimuthal variations of the cosmic ray 
ion acceleration at the blast wave of SN 1006 
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ABSTRACT 

Using radio, X-ray and optical observations, we present evidence for morpho- 
logical changes due to efficient cosmic ray ion acceleration in the structure of 
the southeastern region of the supernova remnant SN 1006. SN 1006 has an ap- 
parent bipolar morphology in both the radio and high-energy X-ray synchrotron 
emission. In the optical, the shock front is clearly traced by a filament of Balmer 
emission in the southeast. This optical emission enables us to trace the loca- 
tion of the blast wave (BW) even in places where the synchrotron emission from 
relativistic electrons is either absent or too weak to detect. The contact disconti- 
nuity (CD) is traced using images in the low-energy X-rays (oxygen band) which 
we argue reveals the distribution of shocked ejecta. We interpret the azimuthal 
variations of the ratio of radii between the BW and CD plus the X-ray and ra- 
dio synchrotron emission at the BW using CR-modified hydrodynamic models. 
We assumed different azimuthal profiles for the injection rate of particles into 
the acceleration process, magnetic field and level of turbulence. We found that 
the observations are consistent with a model in which these parameters are all 
azimuthally varying, being largest in the brightest regions. 
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Introduction 



Collisionless shocks in young supernova remnants (SNRs) are thought to be responsible 
for the production and acceleration of the bulk of Galactic cosmic rays (C Rs) at least up to 



the "knee" (~ 3 PeV) of the CR spectrum (e.g.. iBerezhko fc Volld 120071 ). The theoretical 
mechanism is believed to be the first order Fermi acceleration, also known as diffusive shock 
acceleration (DSA), where electrons, protons and other ions scatter back and forth on mag- 
netic fluctuatio ns or Alfven wave turbulence through the velocity discontinuity associated 
with the shock (IJones &: Ellison! 1 199 It iMalkov &: O'C Druryl I2001L and references therein) . 
The presence of the turbulence is a key ingredient in the process. The higher the turbulence 
level, the higher the scattering rate and the faster the acceleration proceeds, allowing parti- 
cles to reach higher energies. The fast acceleration implies also that the CR pres sure becomes 
signif i cant very qu ickly so that the shock structure is highly nonlinear (e.g 



Ellison et al. 



2000l ; iBaringi 120071 ). The back- reaction of the particles on the shock hydrodynamics in- 
creases the level of magnetic turbulence and the injection rate in turn, i.e., the accelerated 
parti c les create the sca ttering environment by themselves ([Belli Il978l ; iBlandford fc Ostriker 
19781 ; IVolk et al.l 120031 ). Moreover, the production of relativistic particles (extracted from 
the thermal plasma) alon g with their escape from the shock s ystem increases the overall 
density compression ratio (IBerezhko fc Ellison! Il999l ; lBlasill2002l ). 



This theoretical picture leads to several observational consequences. The intensity of 
the synchrotron emission from shock-accelerated electrons will be larger as the acceleration 
becomes more efficient since both the number of relativistic particles and intensity of the 
magnetic field must increase. In the region of very efficient acceleration, the X-ray syn- 
chrotron emission from the highest energy electrons will be concentrated in the form of thin 
sheets behind the blast wave due to synchrotron cooling in the high turbulent magnetic field 
amplifie d at the shock whereas the radio synchrotron emission will be much broader in com- 
parison (lEUison fc Cassam-Chena']]|2005l ). Finally, the morphology of the interaction region 
between the blast wave (BW) and the ejecta interface or contact discontinuity (CD) will 
become much thinner geometrically as the a cceleration becomes more efficient due to th e 
change in the compressibility of the plasma flDecourchelle et al.l |2000| ; lEUison et al.l 120041 ) . 
Because shocks are believed to put far more energy into ions than electrons, this last point, 
if observed, would provide direct evidence for the acceleration of CR ions at the BW. 

The detection of 7-rays with a neutral pion-decay signature resulting from the interac- 
tion of shock-accelerated protons with the ambient matter would provide evidence for CR ion 
acceleration as well. The best chance to see a clear pion-decay signal is when a SNR inter 



acts w ith a dense medium, the prototype case being RX J1713. 7-3946 (ICassam-Chenai et al. 



20041 ). Very high energy 7-rays have been detected from this remnant (lAharonian et al. 
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2007al ) and also from a few others (e.g., lAharonian et al.ll2007bl ; iHoppe et al.l 120071 ). How- 
ever, the pion-decay signal is generally difficult to extract from these observations because 
there are other processes efficiently producing 7-rays but involving CR electrons, namely 
bremsstrahlung and Inverse- Compton scattering off the 2.7 K cosmic microwave background 
and Galactic photons. The nature of the underlying process leading to 7-rays and therefore 



the nature of the emissive particles is sub ject to inten se debate (e.g.. lBerezhko &: Volk!l2006 



Porter et al.ll2006l ; lUchiyama et al.l 120071 ; iPlagal 120081 ) . Besides, in a large number of SNRs 



where bright and geometrically thin X-ray synchrotron-emitting rims are observed (e.g., 
Keper, Tycho, SN 1006), ambient densities are often very low leay ing open the question of 
CR ion production in these SNRs (see also lKatz fc Waxmanll2008l ). 



The idea that the gap betwe en the BW and CD allows quantifying the efficiency of CR 
ion acceleration was proposed by iDecourchellei (120051 ) in the case of Tycho, Kepler and Cas 
A. In Tycho, a remnant where thin X-ray synchrotron-em itting rims are observed all around 
with little variations in brightness, I Warren et al.l (120051 ) showed that the B W and CD are 
so close to each other that they cannot be described by standard hydrodynamical models 
(i.e., with no CR compone nt), hence providing evidenc e for efficient CR ion acceleration in 
this SNR. In a later study, ICassam-Chena'i et al.l (120071 ) showed that the X-ray and to some 
extent the radio properties of the synchrotron emitting rim as well as the gap between the CD 
and B W were consistent with efficient CR ion acceleration modifying the hydrodynamical 
evolution of Tycho. 

Ideally, one would like to find a SNR where both efficient and inefficient particle acceler- 
ation take place so that differential measurements can be made. This offers a great advantage 
for the astrophysical interpretation since a number of uncertainties (e.g., distance and age) 
can be eliminated in the comparison. Furthermore, by directly observing a difference in the 
gap between the BW and CD in regions of efficient and inefficient acceleration, one obtains 
a direct, model-independent, confirmation that CRs do indeed modify the hydrodynamical 
evolution of the BW. SN 1006 provides a valuable laboratory in that regard since the gradual 
variations of the intensity of the radio and X-ray synchrotron emission along the BW point 
to underlying variation of CR ion acceleration. The task of measuring the ratio of radii 
between the BW and CD in the regions of inefficient and efficient CR acceleration will be 
challenging. Nonetheless, SN 1006 is one of the rare objects (perhaps even the unique one) 
where this can be accomplished as we explain below. 

In this paper, we focus our attention on the poorly studied southeastern (SE) quadrant 
of SN 1006, i.e., along the shock front lying between the northeastern and the southwestern 
synchrotron-emitting caps. Our first goal is to measure the ratio of radii between the BW 
and CD there. We can fairly easily follow the BW in the bright synchrotron-emitting caps 
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(presumably the regions of efficient CR acceleration) and the BW can still be traced from 
the Ha emission in the region where both the radio and X-ray synchrotron emissions are 
either too faint to be detected or absent (presumably the regions of inefficient CR acceler- 
ation). As for the CD, it can be traced from the thermal X-ray emission associated with 
the shocked ejecta. Our second objective is to measure the azimuthal variations of the syn- 
chrotron brightness in the radio and several X-ray energy bands. With these observational 
key results in hand, we use CR-modified hydrodynamic models of SNR evolution to constrain 
the azimuthal variations of the acceleration parameters, i.e., the injection rate of particles 
in the acceleration process, strength of the magnetic field, particle diffusion coefficient and 
maximum energy of electrons and ions. 

This paper is organized as followed. In $£2\ we describe the basic characteristics of SN 
1006. In §121 we present the data used in our study. In §H we present the key observational 
results in SN 1006. In §0 we try to relate hydrodynamical models with those results in the 
context of CR-unmodified and -modified shocks. Finally, we discuss our results (§E} and 
present our conclusion (^7j). 



2. Basic characteristics of SN 1006 



SN 1006 was a thermonuclear supernova (SN) widely seen on Earth in the year 1006 
AD. More than a thousand year later, the remnant from this explosion is a huge shell of 30' 
angular size. 

The synchrotron emission is detected from the radio up to the X-rays and dominates in 
two bright limbs - the northeast (NE) and southwest (SW) - where several th in (20" width 
or so ) rims /arcs are running at the periphery, sometimes crossing each other (IBamba et al. 



20031 ; lLong et al.l 120031 ; iRothenflug et al.l 120041 ) . Those apparent ripples are most likely the 



result of the projection of undulating sheets associate d with the shock front. The radio emis- 
sion is well correlated with the nonthermal X-rays ([Reynolds fc Gilmord Il986l ; lLong et al. 
20031 ). Both show an abrupt turn-on with coincident edges in the radial directions, but the 
X-ray emission has more pronounced narrow peaks. This points to a scenario in which the 
rims are limited by the radiative lossejj] in a highly turbulent and ampli fied magnetic field 
faerezhko et al.l I2OO2L I2OO3I : lEllison fc Cassam-Chena] I2OO5I : balletl 120061 1 . Finally, the ob- 
served synchrotron morphology in SN 1006 is best explained if th e bright limbs are polar caps 



with the ambient magnetic field parallel to the shock velocity (jRothenflug et al. 



20041 ). In 



■"■Note however that such a model is incapable of rep roducing the observed shar p turn-on of radio syn- 
chrotron emission at SN 1006's outer edge, as in Tycho (|Cassam-Chena'i et al.ll2007l) . 
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the caps, the maximum energy reached by the accelerat ed particles, as well as their number, 
must be higher than elsewhere (IRothenflug et al.ll2004l ). 



In contrast to the nonthermal X-ray emissi on, the very faint ther mal X-ray emission 
seems to be distributed more or less uniformly (jRothenflug et al.l 120041 ). This is best seen 
in the oxygen band (0.5-0.8 keV). It is, however, difficult to separate the thermal emission 
from the nonthermal emission in the caps. Moreover, it is not clear whether the t hermal 
X-rays should b e attributed to th e shocked ambient medium (jYamaguchi et al.l 120071 ) or the 
shocked ejecta (ILong et al.ll2003l ). but both the over-solar abundances required to fit the 
X-ray spectra in the inner northwest (NW) and NE parts of the SNR and the dumpiness 
on a 30" to 1' scale of the low-energy X-ray emission favor an ejecta origin. The SE region 
seems to be slightly different than the rest of the SNR in terms of ejecta co mposition and 



dum piness too. There is, in p articular, evidence for the presence of cold (IWinkler et al. 



20051 ) and reverse-shock heated (jYamaguchi et al.ll2007l ) iron. I n terms of dynamics, there is 



evidence for ejecta extendin g to/overtaking the BW in the NW (ILong et al.ll2003l ; IVink et al. 
20031 : iRavmond et al.l 120071 ) . 



In the optical, the remnant has a very different morphology. Deep Ha imaging re- 
veals very faint Balmer emission around almost the entire per i phery with a clear arc or 



veals very taint .balmer emission around almost tne entire perrpnery witn a clear arc or 
shock front running from east to south Jwinkler fc Lone] 1 19971 : 1 Winkler et aDl2003h . In- 



terestingly, no synchrotron emission is detected in X-ray in this southeast (SE) region 
whereas extremely faint (an d highly polarized) synchrotron emission is detected in the radio 
([Reynolds fc Gilmordll993l ). Besides, none of the optical emission is of synchrotron origin. 



There is also a clear Balmer 
in gr eat detail in t he optical ( 



2007 ). ultr aviolet jRaymond et al 



ine filament in t 


le NW of SN 1006 which has been observed 


Ghavamian et al. 


2002 




Sollerman et al. 


2003: 


Raymond et al. 



ray bands (ILong et al. 



2003 



1995 



Korreck et al. 



Vink et al. 2003: Acero et al. 



2004 



Hamilton et al.l 120071 ) and X 



20071 ). providing diagnostics for 



the shock speed (V s ~ 2400 — 3000 km/s), ion-electron thermal equilibration at the shock 
(T e /T p < 0.07) and preshock ambient density (0.25 cm -3 < no < 0.4 cm" 3 ). Optical proper 



motions (Winkler et a 



20021 : iHeng fc McCrav 



20031) combined with the shock speed estimate (jGhavamian et al. 
20071 ) along this NW filament leads to a distance to the SNR between 
1.8 a nd 2.3 kpc. Howeve r , because the remna nt is interacting with a denser medium in the 
NW (jMoffett et al.lll993l : iDubner et al.ll2002l ). the shock speed and the ambient d ensity, in 



particular, might be respectively hig her and lower els ewhere (lAcero et al 



2007). For in- 



stance, the broadening of the oxygen (IVink et al.ll2003l ) and silicon (jYamaguchi et al.ll2007l ) 
X-ray emission-lines suggest a shock speed V s > 4000 km/s (although these measurements 
were carried out at different locations in SN 1006). Such high velocities are compatible with 



2001 



recen t proper motion measurements based on radio data (iMoffett et al 
20081 ). Moreover, the lack of observed TeV 7-ray emission (jAharonian et al. 



Revnoso et al. 



20051 ) sets an 
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upper limit on the ambient density of uq < 0.1 cm 3 (IKsenofontov et al.ll2005l ) and the high 
l atitude (~ 50 pc at a distance of 2 kpc) of SN 1006 suggests that n ~ 0.03 — 0.04 cm -3 
rtFerrierdEoOlL 



3. Data 
3.1. X-ray 

We used the most recent Chandra data of SN 1006 (Obsld 3838 and Obsld 4385 up to 
4394) which were obtained in April 2003 with the ACIS-I imaging spectrometer in timed 
exposure and very faint data modes. Eleven pointings were necessary to cover the entire 
extent of the remnant. The final image is shown in Figured! The X-ray analysis was done 
using CIAO software (ver. 3.4). Standard data reduction methods were applied for event 
filtering, flare rejection, gain correction. The final exposure time amounts to ~ 20 ks per 
pointing. 



3.2. Radio 

Radio observations were performed in 2003 (nearly at the same time as the X-ray ob- 
servations) with both the Australia Telescope Compact Array (ATCA) and the Very Large 
Array (VLA). 

The ATCA consists of an array of six 22-m antennae that attain a maximum baseline 
of 6 km in the East- West direction. The ATCA observations occurred on three separate 
occasions in antenna configurations that optimize high spatial resolution: January 24 in its 
6-km B configuration for 12 hours, March 3 in its 6-km A configuration for 12 hours, and 
June 12 in its 750-m C configuration for 7 hours. Observations were made using two 128- 
MHz bands (divided into 32 channels each) centered at 1384 and 1704 MHz. The source 
PKS 1934-638 was used for flux density and bandpass calibration, while PKS 1458-391 was 
used as a phase calibrator. The total integration time on SN 1006 was over 1500 minutes. 

The VLA consists of 27 25-m antennae arranged in a "Y" pattern. Observations with 
the VLA were carried out during 4 hours on January 25th in its hybrid CnD configuration 
(arranged to maximize spatial resolution when observing southern declination sources). Ob- 
servations were made using two 12.5-MHz channels centered at 1370 and 1665 MHz. 3C286 
and 1451-400 were observed for flux and phase calibration. The total VLA integration time 
on SN 1006 was 140 minutes. 
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Data from all of the observations were combined, then uniformly weighted during imag- 
ing to minimize the effects of interferometric sidelobes. The final image, shown in Figure 
121 has a resolution of ~ 6" x 9" and an off-source rms noise of 20 //Jy/beam near the rims. 
Since the angular size of SN 1006 is comparable to the primary beam of both instruments 
(half a degree), corrections were applied to the image to recover the lost flux. After primary 
beam corrections, the total flux density was recovered to within < 5% of the expected value. 
With primary beam correction, noise increases in the faint SE portion. For the study of the 
azimuthal variations of the radio emission at the shock, we use non beam corrected data 



3.3. Optical 



We used the very deep Ha image presented by lWinkler et al.l (120031 ). This image taken 
in June 1998, that is 5 years before the X-ray data, is shown in Figure [3J To compare the 
optical and X-ray images, we had to correct for the remnant's expansion. In the optical, 
proper motions of 0.280 ± 0.008 " yr -1 were measured (from April 1987 to Ju ne 1998) along 



the N W rim where thin nonradiative Balmer-dominated filaments are seen (IWinkler et al. 



20031 ). In the radio, an overall expansio n rate of 0.44 ± 0.13 " yr -1 was measured (from May 
1983 to July 1992. IMoffett et al.lll993 ), although higher values hav e been measured using 
more recent observations ( IMoffett et al.l 120041 ; iReynoso et al.l 120081 ) . This value does not 
include the NW rim where the optical filaments are observed because there is simply little 
or no radio emission there. The observed expansion rate is clearly higher in the radio than 
in the optical and is consistent with the picture in which the BW encounters a localized and 
relatively dense medium in the NW. For simplicity, when using radial profiles, we use the 
value of 0.40" yr -1 , so that the optical profiles are shifted by 2" to be compared with the X- 
ray profiles. As we show below, this correction is negligible compared to other uncertainties 
when trying to locate the fluid discontinuities. 



4. Key observational results 

4.1. Fluid discontinuities 

In this section, we determine the location of the fluid discontinuities - i.e., the BW and 
CD - in the SE quadrant of SN 1006. Our goal is to measure the gap between the BW and 
CD clS db function of azimuth. 

To trace the location of the BW, we use the Ha image of SN 1006 which shows very faint 
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filaments of Balmer emission around much of the periphery of the remnant (Fig. HI top-right 
panel). In the southeastern quadrant, the Ha emission follows a nearly circular arc. We 
extracted radial profiles by azimuthally summing over 4° wide sectors starting from 158° 
(from west) to 302° (in the clockwise direction). The center of these profiles was determined 
to provide the best match to the curvature of the Ha rim. The coordinates of the center are 
(002000; <5j2ooo) = (15 h 02 m 56.8 s , — 41°56'56.6") . This center is very close to the geometrical 
center of the SNR, which, given its overall circularity, may not be too far from the overall 
expansion center which has not yet been determined. In practice we found it very difficult 
to determine the location of the optical rim based purely on a numerical value (e.g., contour 
value or enhancement factor above the local background level) because of the presence of faint 
stellar emission, poorly subtracted stars, and faint diffuse Ha emission across the image. So, 
we identified locations where the Ha emission increased slightly in the radial profiles by eye 
and then further checked (again visually) the corresponding radii on the optical image. We 
associated generous uncertainties with the radii derived from this procedure: 12" outward 
and 24" inward (our process tended to overestimate the radii, hence the asymmetric errors). 

To trace the location of the CD, we use the X-ray image in the low-energy band (0.5 — 0.8 
keV) which contains most of the oxygen lines (Fig. HJ bottom-left panel). While it is not 
clear whether the oxygen emission comes from the shocked ambient medium or the shocked 
ejecta, we consider here an ejecta origin. We discuss and justify this assumption below 
( §5.ip . Another issue comes from the fact that the X-ray emission in the oxygen band may 
contain some nonthermal contribution in the bright limbs. The three-color composite X-ray 
image (Fig. [T| reveals that the oxygen emission (red color) is in fact fairly different from 
the nonthermal X-rays (white color), and appears notably clumpier even in the synchrotron 
limbs (see the eastern limb for instance). This tells us that we can still use the oxygen 
emission to trace the CD in the bright limbs, at least in the azimuthal range we selected 
(158° — 302°) where the synchrotron emission is not completely overwhelming. Radial profiles 
were extracted from a flux image, summing over 1° wide sectors, with the same center and 
azimuthal range as for the optical. In the radial profiles, we selected the radii for which the 
brightness becomes larger than 1.5 x 10 -5 ph/cm 2 /s/arcmin 2 . 

Figure summarizes the above-described measurements and shows the azimuthal varia- 
tions of the radii as determined from the Ha emission (black lines) and low-energy (0.5 — 0.8 
keV, red lines), mid-energy (0.8 — 2.0 keV, green lines) and high-energy (2.0 — 4.5 keV, 
blue lines) X-ray emissions. The cross-hatched lines indicate the location of the bright syn- 
chrotron limbs. They correspond to the places where the contour of the outer high-energy 
X-ray emission (blue lines) can be determined. Between the two synchrotron caps (angles 
between 200° and 270°), there is a clear gap between the BW and CD which is easily seen in 
Figure H] (bottom-left panel). In the bright limbs, however, the BW and CD are apparently 
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globally coincident. Fingers of ejecta (indicated by the stars) are even sometimes visible 
clearly ahead of the Ha emission. In these places, the Ha emission might not be the best 
tracer of the BW location. The mid-energy X-ray emission (green lines) shows that these 
fingers are also present. It is not clear whether this emission traces the BW or, again, the 
ejecta (mostly silicon). The high-energy X-ray emission (blue lines) cannot be used, nor a 
X-ray spectral analysis, to answer this point because of the low number of X-ray counts in 
this region. 

Figure [6] shows the ratio of radii for the Ha and low-energy X-ray emission, which to 
first approximation corresponds to the ratio of radii between the BW and CD, Rbw/Rcd- 
With increasing azimuthal angles the ratio of radii increases from values near unity (within 
the northeastern synchrotron-emitting cap) to a maximum of Rbw/Rcd — l-10too4 before 
falling again to values near unity (in the southwestern synchrotron-emitting cap). Over the 
entire azimuthal region where the synchrotron emission is faint, the azimuthally averaged 
ratio of radii is Rbw/Rcd — 1-04 ± 0.03. In the regions within the synchrotron rims, 
Rbw/Rcd ^ LOO. 



4.2. Synchrotron emission at the blast wave 

In this section, we measure the azimuthal variations of the synchrotron flux extracted 
at the BW and we investigate how these variations compare when measured at different 
frequencies. 

For that purpose, we use radio and X-ray data (Fig. m left panels). We extracted the flux 
in a 30"-wide region behind the BW along the SE rim. To define the BW location, we used the 
well-defined radii obtained from the high-energy X-ray image for the northeastern (9 < 200°) 
and southwestern (9 > 270°) regions while, in the remaining SE quadrant (200° < 9 < 270°), 
we used the radii determined from the Ha image since no synchrotron emission is detected 
there. 

In Figured (top panel), we show the azimuthal variations of the projected brightness in 
the radio (1.5 GHz, black lines), and several X-ray energy bands: 0.5 — 0.8 keV (red lines), 
0.8 — 2.0 keV (green lines) and 2.0 — 4.5 keV (blue lines). The radio brightness was multiplied 
by 2 x 10~ 5 to make it appear on the same plot with the X-ray data. When moving from the 
northeastern or southwestern synchrotron-emitting caps toward the faint SE region, we see 
that the projected brightness gradually decreases in both radio and X-ray bands. The radio, 
medium- (green lines) and high-energy (blue lines) X-rays show mostly the variations of the 
synchrotron emission, while the low-energy X-rays (red lines) show mostly the variations of 



-10- 



the thermal (oxygen) emission, but may contain some nonthermal contribution as well. 

In Figure [7] (bottom panel), we show the same azimuthal variations but rescaled roughly 
to the same level in the synchrotron-emitting caps. The logarithmic scale shows that the 
brightness contrast between the bright caps and faint SE region is of order 20 to 50 in the 
radio and X-rays. In fact, this value might even be larger since the radio and high-energy 
X-ray fluxes in the faint SE quadrant are compatible with no emission. We note that the 
azimuthal variations of the radio synchrotron emission are much less pronounced than those 
in the X-rays. 



5. Relating hydrodynamical models to the observations 

In the previous section, we have found evidence for azimuthal variations of the separation 
between the BW and CD and synchrotron brightness just behind the BW (©. Both the 
separation and the brightness are correlated: the brighter the BW, the smaller the separation. 
We measured in particular an unexpectedly small separation between the BW and CD (~ 
1.04) in the SE quadrant of SN 1006 where little to no synchrotron emission is detected. 
These measurements depend on the assumption that the oxygen emission comes from the 
shocked ejecta. After justifying this approach below ( §5.11) . we first investigate the different 
scenarios that could potentially lead to a small ratio of radii between the BW and CD 
assuming no CR acceleration at the shock ( §5.2)) . Finally, we interpret the observed azimuthal 
variations of the ratio of radii and synchrotron brightness using CR-modified hydrodynamic 
models ( §5.3)) . 



5.1. Does the oxygen come from the ejecta? 

A strong assumption made in the measurement of the ratio of radii between the BW 
and CD in SN 1006 is that the thermal X-ray emission (from the oxygen) is associated with 
the shocked ejecta. 

There is, however, still the possibility that the oxygen emission is associated with the 



shocked ambient medium. lYamaguchi et all (120071 ) have shown that the integrated X-ray 



spectrum of the SE quadrant can be described by a combination of three thermal plasmas 
in non-equilibrium ionization and one power-law component. One of the thermal compo- 
nents, assumed to have solar abundances and therefore associated with the shocked ambient 
medium, was able to produce most of the observed low-energy X-rays (in particular the Ka 
lines from O VII, O VIII and Ne ix). The two other thermal components, with non-solar 



- 11 - 



abundances, were able to reproduce most of the Ka lines of Si, S and Fe, and were attributed 
to the shocked ejecta. 

There are several arguments from an examination of the high resolution Chandra images, 
however, that favor an ejecta origin for the low-energy X-ray emission. The three-color 
composite X-ray image (Fig. [T|) shows that the spectral character of the X-ray emission in the 
SE quadrant (where there is little to no synchrotron emission) does not c hange appreciably 



as a f unction of radius behind the BW, which would be expected if the lYamaguchi et al. 



(120071 ) picture were correct. The lack of virtually any X-ray emission between the BW and 
CD c an be explain ed by the expected low density of the ambient medium (~ 0.03 — 0.04 



cm 



Ferrierdl200ll ) which can reduce both the overall intensity and the level of line emission 
(due to strong nonequilibrium ionization effects). If the soft X-ray emission were from the 
shocked ambient medium its brightness should in principle rise gradually behind th e BW as 
a result of the projection of a thick shell onto the line-of-sight (IWarren et al.ll2005l ). but the 
radial profile of the low-energy X-ray emission (not shown) in fact turns on rather quickly 
at the periphery. 

Another argument in favor of an ejecta origin is the observed dumpiness of the oxygen 
emission throughout the SE region. There is also the protuberances seen right at the edge 
of SN 1006, which are suggestive of Rayleigh- Taylor hydrodynamical instabilities that are 
expected at the CD (see Fig. [1]). Note that the outer edge of the long, thin filament of 
X-ray emission in the NW - which overlaps the brightest Ha emission and therefore is likely 
to be at least partly due to shocked ambient medium - is much smoother than the edge 
of the remnant in the SE. Finally, high emission measures of oxygen in the shocked ejecta 
at young dynamical SNR ages are compatible with most thermonuclear explosion models. 
Essentially, all 1-D and 3-D deflagration models, and all delayed detonation and pulsating 
delayed detonation models have oxygen in the outer layers (but not prompt detonations or 
sub-Chandrasekhar explosions) (IBadenes et al.ll2003l ). 



5.2. The small gap in the faint SE region 

We have presented evidence in §4.11 for a small ratio of radii Rbw/Rcd ~ 1-04 between 
the BW and CD in the SE quadrant of SN 1006. This is also the location where the radio and 
X-ray nonthermal emission are either very weak or undetected. The most straightforward 
explanation for this lack of synchrotron emission is an absence of efficient CR acceleration at 
the BW here. In the following sections, we first investigate whether pure 1-D hydrodynamical 
models for SNR evolution (i.e., without CR acceleration) can reproduce such a small ratio 
of radii. We consider the role of ambient density, ejecta profile, and explosion energy on the 



-12- 



size of the gap between CD and BW. Then we consider the effects of 3-D projection on the 
hydrodynamically unstable CD. 



5.2.1. Predictions from standard hydrodynamical models 



To understand how a small gap can be produced in a region where there is apparently no 
evidence for efficient particle acceleration, we used standard one dimensional (1-D) spher- 
ically symmetric numerical hydrodynamical simulations that follow the interaction of the 
ejecta with the ambient medium. These simulations do not include any CR component. 

Key inputs to these simulations are the initial density profiles of the ejecta and ambi- 
ent medium. In the following, we consider a uniform density in the ambient medium and 
two different initial ejecta density profiles: exponential and power-law. The exponential 
form has been shown to be m ost representative of explosion models for thermonuclear SNe 
( IDwarkadas &: Chevalierlll998l ). but we also present results with a power-law distribution to 
quantify the impact of the shape or compactness of the ejecta profile on Rbw/Rcd- 

In general, an exponential profile is expected for thermonuclear SNe because the explo- 
sion is driven by the continuous release of energy from the burning front as it propagates 
through the star, while the power-law profile (more compact) is expected for core collapse 
SNe because the explosion is driven by a ce ntral engine (core bounce ), and the shock loses 
energy as it propagates through the star (IMatzner &: McKed Il999l ). Power law profiles 
with n = 7 have been us ed to represent 1-D deflagration models (in particular, model W7 
from lNomoto et al.lll984T), while exponential profil es are more adequate to represent delayed 
detonation models (IDwarkadas fc Chevalier! Il998l ). When comparing these two analytical 
profiles, it is worth noting that they transmit momentum to the ISM in a different manner. 
The more compact power law profile is a more efficient piston for ISM acceleration, and will 



lead to smaller gaps than the less compact exponential profile (IDwarkadas fc Chevalierlll998 
Badenes et aDl2003h . 



The hydrodynamical simulations provide the radii of the BW and CD at any time for a 
given ejecta profile, ambient medium density, and explosion energy. In Figure [BJ we quantify 
the impact of these three contributions on the ratio of radii Rbw/Rcd- In particular, we 
plot the ratio obtained for a wide range of ambient medium density (0.001 cm" 3 < n < 
30 cm -3 ) using 1-D exponential (EXP) and powerlaw (PL7) ejecta profiles with different 
kinetic energies. The cross-hatched domains define the range of ratio of radi@ observed in 



2 In Figure [8l it is important to keep in mind that the maximum value of the ratio found in Tycho comes 
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SN 1006 and Tycho (Rbw/ Rcd < 1.10) and range of typical ambient densities (0.01 cm -3 < 
n < 0.06 cm -3 in SN 1006; 0.1 cm" 3 < nn < 0.6 cm~ 3 in Tycho) consistent wit h constraints 
derived from the observations (lAcero et al.l 120071 ; ICassam-Chenai et al.l 120071 ) . It is clear 
that the 1-D hydrodynamical simulations are unable to reproduce a ratio of radii as small 
as the one observed in both SNRs. We note however that such a comparison between the 
observations and the models is not straightforward since our measurements are likely to be 
affected by projection and other effects as we detail next. 



5.2.2. Three dimensional projection effects 



One of the limitations of our numerical simulations is the one dimensionality. Hydrody- 
namical simulations in 2-D or 3-D show that hydrodynamical (Rayleigh-Taylor) instabilities 
project pieces of ejecta ahead of the 1-D CD. In these simulations, t he outermost pieces of 



ejecta reach half way of the gap between the 1-D CD and BW (e.g., IChevalier et al.lll992 



Wang &: Chevalier! l200ll ). Therefore, because of such protrusions, the line-of-sight projected 
CD radius, -Rcd, will appear larger than the true average CD radius, -Rcd- On the other 
hand, because the BW is not as highly structured as the CD interface (it is not subject to 
Rayleigh-Taylor instabilities), it is reasonable to assume that the true average BW radius, 
-Rbw, is ver y close to the projected value, -Rbw- It follows that the ratio of projected radii 
between the BW and CD, -R B w/-Rcd, will be smaller than the true average ratio, -Rbw /Rcd- 
Quantifying the effects due to projection is of considerable importance for interpreting the 
ratio of radii between the BW and CD in the context of particle acceleration^. 

There are several attempts aimed at estimating the projection correcting factor, £, 
where £ is defined as Rqb = (1 + -Rcd- Based on the projection with ejec ta protrusions 



determ ined from a power-spectrum analysis done at the observed CD in Tycho. I Warren et al. 



(120051 ) estimated an amount of bias of £ 



proj 



6%. Based on the projection of a shell of 



shocke d eje cta with protrusion s calcu lated in a 2-D hydrodynamical simulations, iDwarkadas 
(120001 ) and Wang fc Chevalier! (1200 ll ) found a slightly larger value of £ pro j ~ 10%. Taking 
£ = 10%, the predicted ratio of projected radii, -Rbw /Rcd = -Rbw /Rcd / (1 + £), becomes 
1.18/(1 + 0.1) ~ 1.07, where 1.18 is the ratio of radii found at an age of 1000 years using 



from a presumably efficient particle acceleration shock region while in SN 1006, it comes from an inefficient 
particle acceleration region. 

3 It is not even clear how the 1-D CR-modified hydrodynamical results on the ratio of radii need to be 
adjusted given mu l ti-D R ayleigh-Taylor instability effects. Preliminary studies on this have been done by 



Blondin fc Ellisonl (|200ll ) but this work does not fully quantify all relevant effects. 
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1-D power-law ejecta profiles assuming a kinetic energy of the explosion of 10 51 ergs and an 
ambient medium density of 0.03cm -3 . Using exponential ejecta profiles, we find -Rbw/-Rcd = 
1.14 because they produce higher ratio of radii (i.e., -Rbw/-Rcd = 1.25). In either case (power- 
law or exponential ejecta profiles), the average ratio (resp. 1.07 or 1.14) is still larger than 
the average value of 1.04 measured in the SE quadrant of SN 1006 (cf. §4. ID . 

There are several possibilities to explain such discrepancy between the models and the 
observations. For instance, the value of £ pro j could be larger in 3-D than in 2-D. Indeed, 
models show that hydrodynamical in stabilities can gro w considerably faster (by ~ 30%) and 
penetrate further in 3-D than in 2-D (IKane et al.ll2000l ). There can also be a certain amount 
of inhomogeneity in the ejecta density distribution. Spectropolari metric observat i ons o f 
thermon uclear SNe do show that ejecta are clumpy on large scales ([Leonard et al.l 120051 ). 
Besides, IWang fc Chevalier! (120011 ) have shown that high density clump traveling with high 
speed in the diffuse SN ejecta can reach and perturb the BW. This requires a high density 
contrast of a factor 100. The density contrast is certainly not as high in SN 1006 because 
otherwise we would see it in the low-energy X-ray and/or in the optical image, but may 
be non-negligible. According to the previous calculation, the ejecta clumping is required to 
contribute an additional £ c i ump ~ 3 — 10% factor (so that £ = £ pro j +£ c ium P ) in order to make 
the observed and predicted ratios of radii agree. Finally, although we could not come up 
with a reasonable solution, we cannot exclude the possibility of a complex 3-D geometry in 
the SE quadrant of SN 1006 where the outermost extent of the shocked ejecta would not 
be directly connected to the shock front emission in the hydrodynamical sense (i.e., the CD 
and BW would corresponds to different parts of the remnant seen in projection). 



5.3. Azimuthal variations of the CR ion acceleration 

In this section, we try to explain both the observed azimuthal variations of the ratio of 
radii between the BW and CD (gTQ and synchrotron flux at the BW (g^D in SN 1006 in 
the context of CR-modified shock hydrodynamics where the efficiency of CR acceleration is 
varying along the BW. 

In the following, we use CR-modified hydrodynamic models of SNR evolution which 
allow us for a given set of acceleration parameters to calculate the ratio of radii, Rbw/R-cd, 
and the synchrotron flux, F u , at any frequency v. Relevant parameters are the injection 
rate of particles into the DSA process, the magnetic field strength and the particle diffusion 
coefficient (or equivalently the turbulence level). Our goal will not consist in finding the exact 
variations of these parameters along the shock by fitting the data, but rather to provide a 
qualitative description of the variations based on the theory to verify whether the model 
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predictions are consistent with the observations or not. 



5.3.1. Conceptual basis for modeling 



Here, we explicate the various assumptions concerning the spatial distribution of the 
ambient density and magnetic field before running our CR-modified hydrodynamic models. 
The knowledge of the remnant's environment is important since it affects its evolution and 
also its emission characteristics. 

First, we assume that the ambient density is uniform. This is suggested by the quasi- 
circularity of the Ha filament in the SE of SN 1006 and a reasonable assumption for remnants 
of thermonuclear explosion (IBadenes et al.l 120071 ). Second, we assume that the ambient 
magnetic field direction is oriented along a preferred axis. In SN 1006, this would correspond 
to southwest-northeast axis which is the direction parallel to the Galactic plane (following 
Rothenflug et al.ll2004j ). The magnetic field lines are then parallel to the shock velocity in 
the brightest part of the synchrotron caps in SN 1006 and perpendicular at the equator. 
The angle between the ambient magnetic field lines and the shock velocity will determine 
the number of particles injected in the PSA process, i. e ., the injection ra te. The smaller the 
angle, the larger the injection rate (lEllison et al.lll995t IVolk et al.ll2003l ). 



With this picture of an axisymmetry around the magnetic field orientation axis in mind, 
the azimuthal conditions at the BW are essentially spatially static, the evolution is self- 
similar and hence temporal evolution of gas parcels can be followed with a 1-D spherically 
symmetric code. We will run such a code with specific initial conditions at each azimuthal 
angle along the BW (viewing the SNR in a plane containing the revolution axis of the 
magnetic field). We implicitly assume a radial flow approximation in our approach, i.e., that 
each azimuthal zone is sufficiently isolated from the others so that they evolve independently. 

We run the self-similar models assuming a power-law density profile in the ejecta with 
an index of n — 7, an ejected mass and kinetic energy of the ejecta of 1.4 M and 10 51 
ergs, respectively, which are standard values for thermonuclear SNe. We assume that the 
SNR evolves into an interstellar medium whose density is uq = 0.03 cm -3 and pressure 2300 
K cm -3 . For such a low density, the use of the self-similarity is well justified as shown in 
Figure [HI the ratio of radii Rbw/Rcd obtained with a 1-D purely numerical hydrodynamic 
simulation with the same initial conditions remains indeed more or less constant even until 
the age of 1000 yrs (solid yellow lines). Other input parameters will be specified later. First, 
we describe the CR acceleration model used with the hydrodynamical model. 
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5.3.2. CR acceleration model 



The models consist of a self-similar hydrodynamical calculation coupled with a nonlinear 
diffusive shock acceleration model, so that the back-r eaction of the accelerated particles at the 
BW is taken into account fjDecourchelle et al.ll2000l ). For a given injection rate of protons, 
r/inj, far upstream magnetic field, B u , and diffusion coefficient, D, the shock acceleration 
model determines t he shock jump conditions and the particle spectrum from thermal to 
relativistic energies (IBerezhko fc Ellison! Il999l ). This allows us to compute the ratio of radii 
and, by following the downstream evolution of the particle spectrum associated with e ach 
fluid element, the synchrotron emission within the remnant (jCassam-Chena'i et al.ll2005l ). 



The CR proton spectrum at the BW is a piecewise power-law model with an exponential 
cutoff at high energies: 

f p (E) = a E~ v ^ exp (-£/£ p , max ) , (1) 

where a is the normalization, V is the power-law index which depends on the energy E, 
and -E^max is the maximum energy reached by the protons. Typically three distinct energy 



regimes with different T values are assumed (IBerezhko fc Ellison! Il999l ). The normalization, 
a, is given by: 

a = (2) 
47r Pinj 

where n in j = ?7i n j n (r tot /r sub ) is the number density of gas particles injected in the acceler- 
ation process, with r to t and r su b being the overall density and subshock compression ratios, 
<7sub = 3 r su b/(r su b — 1) and p in j is the injection momentum. Here, p- m j = X m p c S 2 where A is 
a parameter (here fixed to a value of 4) which encodes all the complex microphysics of the 
shock and c S 2 is the s ound speed in the shock-h eated gas immediately downstream from the 
subshock (see §2.2 in IBerezhko fc Ellison! Il999l ). We note that because there is a nonlinear 
reaction on t he system due to the injection, the parameter ?7 in j should be in fact related to 
Pinj (see §5 in lBlasi et al.ll2005l . for more details). 



The CR electron spectrum at the BW is determined by assuming a certain electron- 
to-proton density ratio at relativistic energies, K cp , which is defined as the ratio between 
the electron and proton distributions at a regime in energy whe re the protons are already 
relativistic but the electrons have not yet cooled radiatively (e.g.. Ellison et al.ll2000l ). In the 
appropriate energy range, the CR electron spectrum is then: 



fe(E) 



a K cp E~ T ^ 



exp (-E/E e 



(3) 



where -E c ,max is the maximum energy reached by the electrons, which could be eventually 
lower than that of protons (E Ptmax ) through synchrotron cooling of electrons depending on 
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the strength of the post-shock magnetic field. K ep is left as a free parameter. We will obtain 
typically K ep ~ 10" 4 - 10" 2 . 



The maximum energies E pmax and E e>max contain information on the limits of the ac- 
celeration. They are set by matching either the acceleration time to the shock age or to the 
characteristic time for synchrotron losses, or by matching the upstream diffusive length to 
some fraction, £ s , of the shock radius, whichever gives the lowest value. We set C, s = 0.05, 
a value that allows to mimic the effect of an expandin g spherical shock as compared to th e 
plane-parallel shock approximation assumed here (see iBerezhkol Il996l ; Ellison et al.ll2000l ). 
A fundamental parameter in determining the maximum energy of particles is the diffusion 
coefficient, D, which contains i nformation on the level of turbulence and encodes the scat- 
tering law fjParizot et al.ll2006l ). For the sake of simplicity, we assumed the Bohm regime 
for all particles (i.e., diffusion coefficient proportional to energy). We allow deviations from 
the Bohm limit (which corresponds to a mean free path of the charged particles equal to 
the Larmor radius, which is thought to be the lowest possible value for isotropic turbulence) 
via the parameter ko defined as the ratio between the diffusion coefficient, D, and its Bohm 
value, Dq. Hence ko = D/Dq is always > 1 and ko — 1 corresponds to the highest level 
of turbulence. The smaller the diffusion coefficient (or ko), the higher the maximum energy 
(the maximum energy scales as 1/vfco i n the radiative loss case and as 1/fco m the age- and 
escape-limited cases). 



5.3.3. Heuristic models 

In order to make predictions for the azimuthal variations of the ratio of radii and syn- 
chrotron flux, we must make some assumptions about the azimuthal variations of the input 
parameters used in our CR-modified hydrodynamical models. Relevant input parameters 
are the injection rate (rj- m j), the upstream magnetic field (B u ), the electron-to-proton ratio 
at relativistic energies (K ep ) and the level of turbulence or diffusion coefficient relative to 
the Bohm limit (ko). 

To understand the influence of these four basic input parameters, we begin first with 
a heuristic discussion to set the stage for later detailed and more physical models ( §5.3.5[) . 
In Figure [9j we show three models where we vary the two most important parameters, 77^ 
and B u , with azimuthal angle (-Kep is held constant and ko is set to 1 at all angles). Model 
1 (left panels) has 77^ varying from 10~ 5 to 10 -3 , while B u is held constant at the value 
25 /iG. Model 2 (middle panels) has constant r] m j = 10 -4 and B u varying from 3 /iG to 
25 /iG. Model 3 (right panels) is a combination of models 1 and 2 and has 77i n j varying from 
10~ 5 to 10~ 3 and B n varying from 3 /iG to 25 /iG. 
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Note that in all these models, the functional form of the profiles when r)- m j and B u vary 
is rather arbitrary and just used for illustration purposes (for completeness we note that the 
angular range shown corresponds to the data extraction region in Fig. [7]). Moreover, the 
upstream magnetic field, B u (second panels, solid lines), is allowed to vary with azimuthal 
angle because we implicitly assume that it can be significantly amplified by the CR-streaming 
instability. Our model does not include self-consistently the magnetic field amplification 
believed to occur at SNR shocks, but is provided with a simple compression. Assuming that 
the magnetic turbulence is isotropic ahead of the shock, the immediate post-shock magnetic 
field, Bd (dashed lines), will then be larger than upstream by a factor r# = B&/ B u = 
a/(1 + 2 rjf ot ) /3. This relation is assumed at each azimuthal angkfl 

We show the predicted profiles of the ratio of radii between the BW and CD (third pan- 
els) and synchrotron brightness projected along the line-of-sight in the radio and different 
X-ray energy bands (bottom panels). We explain later how those brightness profiles were 
precisely constructed ( §5.3.5[) . Clearly, model 2 does not predict the expected correlation 
between the ratio of radii and the synchrotron brightness and hence can be rejected imme- 
diately. On the other hand, models 1 and 3 lead to the expected correlation: the smaller 
the BW/CD ratio, the brighter the synchrotron emission. Comparing models 1 or 3 with 
model 2 allows us to see that only an azimuthal variation of the injection rate leads to the 
appropriate variation in the ratio of BW to CD radii (third panels). Comparing model 1 
with model 3 shows the impact of the magnetic field variations through synchrotron losses 
on the azimuthal profiles of the X-ray synchrotron emission (bottom panels). 

Although models 1 and 3 seem to describe the key observational constraints fairly well, 
they suffer weaknesses as regards their underlying astrophysical premises. In model 1, it 
is difficult to justify why the magnetic field should be enhanced in the region of inefficient 
particle acceleration (77^ < 6 x 10~ 5 ). The model proposes a post-shock magnetic field 
there of ~ 80 fiG, when a value of 10 — 15 /xG would more likely reflect the value of the 
compressed ambient field. In model 3, it is difficult to understand why the turbulence level 
would saturate (i.e., ko — 1) everywhere. Hence, we are led to develop more elaborated 



4 Here, the compressi on does not dep end on the angle betwee n the upstream magnetic field and shock 
velocity as considered in iRevnolds ( 1998 ) or Orlando et al. (2007). Instead we assume that after the ampli- 



fication process, the magnetic field becomes entirely turbulent, so its originally ordered character has been 
lost. 

5 Looking at the BW to CD ratio in model 2, it naively seems surprising to have a less modified shock 
(i.e., higher BW to CD ratio) where the magnetic field is larger. This effect is due to the presence of Alfven 
heating in the precursor (see Berezhko fc Ellison 1999h . although th ere is debate about w hether this is the 
most efficient mechanism for turbulent heating in the precursor (see lAmato fc Blasill2006l ). 
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models where and/or K ep are allowed to vary with azimuthal angle. A variation of either 
fco or A' ep will not significantly affect the ratio of BW to CD radii, but will only modify 
the profiles of the synchrotron emission. In fact, a model with varying K ep can be already 
excluded because this quantity changes the radio and X-ray synchrotron intensities in the 
same way while the observed radio and X-ray profiles (Fig. [7]) vary with azimuth in different 
ways. Varying ko (as we show in §5.3.5p does allow us to modify the relative azimuthal 
profiles. 

5.3.4- Can CR-modified models provide Abw/Acd = 1.00? 

Although CR-modified hydrodynamical simulations predict a ratio of radii -Rbw/-Rcd 
smaller than standard hydrodynamical simulations, they still predict a lower limit on that 
ratio that is strictly larger than unity. For instance, Figure [9] (third-right panel) shows that 
the modeled ratio of radii where the synchrotron emission is strong is ~ 1.13. This value 
was obtained assuming an injection rate of 10~ 3 and an ambient magnetic field of 25 /xG. 
Keeping all the same inputs, but reducing the magnetic field value to 3 yuG, we obtain the 
lowest possible value of Rbw/Rcd = 1.06. 

On the other hand, the observations of SN 1006 show that the ratio of radii is of order 
1.00 in the synchrotron-emitting caps, a value that is not possible to obtain from our CR- 
modified model under the DSA framework. However, as already mentioned before in our 
discussions of the small gap between the BW and CD in the region of inefficient acceleration 
( §5.2p . we need to consider projection effects due to hydrodynamical instabilities at the CD 
which can reduce the gap by 6 — 10%. This does not fully account for the difference with the 
model and, therefore, we may need to invoke another effect, such as clumping of the ejecta, 
to fully explain the offset between the observation and models. It is important to note that 
in both the regions of efficient and inefficient CR acceleration, we require essentially the same 
numerical factor to bring the modeled ratio of radii (for the appropriate model in each case) 
into agreement with the observed ratio. In the following, we introduce an ad-hoc rescaling 
of the ratio of radii by decreasing the modeled values by a constant factor (~ 13%) at all 
azimuthal angles. By doing this, we implicitly assume that this factor contains the total 
contribution of effects that are unrelated to the CR acceleration process. The numerical 
value used depends, of course, on the particular SN explosion model used. For example, the 
factor would need to be increased slightly in the case of an exponential ejecta profile (cf., 

Fig. ED. 
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5.3.5. Toward a good astrophysical model 

In Figure [101 we present in more detail what we believe is a reasonable astrophysical 
model. This model assumes that ?7i n j increases from a value of 5 x 10~ 5 to a value of 10~ 3 while 
at the same time B u increases from a value of 3 fiG to a value of 20 fiG and ko increases from 
an arbitrary value of 100 to 1 (K ep is kept constant at all angles). The azimuthal variations 
of these parameters is consistent with the picture of DSA where, once the magnetic field 
fluctuations at the shock exceed the background ISM fluctuations, any low injection rate 
leads to growth of the turbulent m agnetic field in th e upstream region, which in turns leads 



to an increase of the injection rate (jVolk et al.ll2003l ). The acceleration efficiency, defined as 



the fra ction of total energy flux crossing the shock that goes into relativistic particles (see 



3.3 in iBerezhko fc Ellisonl Il999f ). is 51% in the case with r)- m j = 10 3 and 7% in the case 



with 77inj = 5 x 10 5 . 

Note that the form of the azimuthal profiles of ?7i n j and B u is again arbitrary (here a 
squared-sine variation) and just used for illustration purpose. Our goal is not to find the exact 
form of the profiles based on an accurate fit of the model predictions to the observations, 
but rather to describe qualitatively how the acceleration parameters vary along the BW in 
SN 1006. The extremum values for r] m j were chosen so that the predicted variations of the 
ratio of radii Rbw I -Rcd roughly matches the observations. A maximum (minimum) value of 
Tfrnj between 10~ 4 and 10~ 2 (resp. 10~ 5 and 5 x 10~ 5 ) will not significantly affect our results 
concerning the profile of Rbw/Rcd, but differences will be seen in the absolute intensity 
in the synchrotron emission. The maximum value of B u was adjusted so that it yields an 
immediate postshock value, B^ ~ 80 fiG (panels e), consistent with the value derived from 
the thickness of the X-ray synchrotron-emitting rims in SN 1006, under the assumptio n that 



this t hickness is limited by the synchrotron losses of the highest energy electrons (e.g., [Ballet 



20061 ). Finally, we adjusted the minimum value of B u in the faint region so that the predicted 
and observed radio synchrotron flux are roughly the same; here we show profiles with a value 
of B u = 3 fiG, i.e., the typical ISM value. 

In Figure [TU1 (panel g), we show the predicted azimuthal profile of ratio of radii between 
the BW and CD (dashed lines). The ratio of radii was normalized (solid lines) so that it 
roughly equals 1.00 in the region of efficient CR acceleration as in the observations. This 
normalization procedure (which we justified in §5.3.4p works remarkably well. The variations 
of the ratio of radii reflects nothing but the variations of the overall density compression ratio, 
ftot (panel d). For low injection rates (77^ < 6 x 10~ 5 ), the compression ratio comes close to 
the value obtained in the test-particle case (i.e., r tot = 4) and as the injection rate increases, 
the plasma becomes more compressible with r tot ~ 6. However, once the magnetic field has 
been considerably amplified (B u > 10 fiG), r tot starts to slightly decrease due to the heating 
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of the gas by the Alfven waves in the precursor region (see iBerezhko fc Ellison! Il999l ) 



In Figure [TU] (panel h), we show the azimuthal profiles of the synchrotron emission in the 
radio and several X-ray energy bands. To obtain these profiles, we first computed the radial 
profile of the synchrotron flux (calculated in a given energy band) at each azimuthal angle. 
Then, we projected the flux along the line-of-sight and finally extracted the projected flux in 
30"-wide regions behind the BW (assuming a distance of 2 kpc) as we did in the observations 
of SN 1006 ( §4.2p . In this procedure, we assumed a spherically symmetric distribution for 
the emissivity. This is a reasonable assumption as long as we restrict the analysis close to 
the BW, considering that the bright limbs in SN 1006 are in fact polar caps, i.e., where the 
emissivity distribution is axisymmetric. Finally, we fixed the radio flux in the synchrotron 
caps to be the same as in the observations by adjusting the electron-to-proton density ratio 
at relativistic energies, K ep . This is obtained for a value of K ep = 7 x 10~ 4 . 

Relaxing the assumption of the Bohm diffusion (i.e., ko) at each azimuthal angle sig- 
nificantly modifies the azimuthal profiles of the synchrotron flux in the X-ray band (while 
the radio profile is unchanged) compared to the model 3 of Figure M in which & was always 
equal to 1. This is because an increase of k decreases the maximum energy of electrons, 
-Ee.max- Since the synchrotron emission is very sensitive to the position of the high-energy 
cutoff in the electron distribution (i.e., -E^max), all profiles do not show anymore a plateau 
in the region of efficient CR acceleration but are now gradually decreasing (panel h). Be- 
sides, the synchrotron brightness in the X-ray bands now decreases faster than in the radio 
(panel i). This is roughly consistent with what is observed in SN 1006. We note that in this 
model, the maximum energy of both protons and electrons (panel f) decrease rather quickly 
as we move toward regions of inefficient CR acceleration as opposed to the previous model 3 
with k = 1 where the maximum energy of electrons remained approximat ely constant (not 



shown). Overall, our results from this model are consistent with those of iRothenflug et al. 



(120041 ) which were based on a different approach (i.e., measure of the azimuthal variations 



of the cutoff-frequency all around SN 1006) using XMM-Newton observations. 



5.3.6. Other azimuthal profiles for rj- m j, B u and k 

We have found that a simple model in which the injection rate of particles in the acceler- 
ation process (r/mj), amplified magnetic field (B u ) and level of turbulence (k ) were gradually 
increasing around the shock front was able to provide predictions in a good agreement with 
the observations for the azimuthal variations of the ratio of radii and radio and X-ray syn- 
chrotron fluxes. The input parameters and their azimuthal variations were chosen based 
on the physical picture where the injection of particles and amplification of the magnetic 
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turbulence are coupled. 

One can imagine however more complicated scenario and hence more complicated input 
profiles for the acceleration parameters (if m j, B u , k ). For instance, the coupling between 
the injection and the turbulence which occurs at the beginning of the acceleration process 
may end at some point when the amplitude of the magnetic field fluct uations become s o high 



that their further growth is prevented by strong dissipation processes ( iVolk et al.ll2003l ). The 
injection, magnetic field and level of turbulence may in fact saturate (over some azimuthal 
range) as we reach regions of very efficient CR acceleration. Another example of scenario can 
be imagined if the shock velocity is not uniformly distributed along the shock (as opposed 
to what we have assumed in our models) due to a asymmetrical explosion. Such a scenario 
could explain why the bright synchrotron-emitting caps in SN 1006 seem to have larger radii 
than the faint regions. If so, the injection could be larger in the region of higher velocity 
while the level of turbulence and magnetic field could potentially saturate there. In general, 
these new scenarios will not lead to a significant modification of the profile of the ratio 
of radii. However, we do expect some changes in the azimuthal profiles of the radio and 
X-ray synchrotron flux. They will generally be flatter in the regions of very efficient CR 
acceleration. 



6. Discussion 

6.1. Polar cap morphology 

The model in which the injection rate of particles, magnetic field and level of tur- 
bulence are all varying as a function of azimuth was successful in reproducing the overall 
azimuthal variations of the radio and X-ray synchrotron flux observed along the BW of SN 
1006 (Fig. [TUT) . In fact, it is actually possible to build a 2-D projected map of the radio and 
X-ray synchrotron morphology comparable to the radio and X-ray images of SN 1006 (Figs. 
[Hand [2]). This requires us to know the t hree-dimensional distrib ution of the emissivity. Here, 



we consider polar caps as suggested by iRothenflug et al.l (120041 ) . 



To simplify the calculation, we consider that the radial profiles of the emissivity obtained 
from the different values of injection, magnetic field and turbulence level have an exponential 
form characterized by a maximum emissivity at the shock and width AR over which the 
emissivity decreases (see Appendix |A]) . We also fixed the BW radius, Rbw, to be the same 
at each azimuthal angle where the injection and magnetic field vary, although it slightly 
depends on the CR acceleration efficiency. These are good approximations. 



In the case of the radio synchrotron emission, the width AR does not depend on the 
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azimuthal angle (from the pole to the equator), and we have typically AR = AR/Rbw — 
0.01. The projected radio morphology is shown in Figure [TT] (top-right panel). It has a 
bipolar limb-brightened morphology as is observed in SN 1006. 

In the case of the X-ray synchrotron emission, the width AR is expected first to increase 
(starting from the equator) as particles reach higher energies and the emission builds up, and 
then to slightly decrease when the effect of radiative losses starts to become important due 
to the larger post-shock magnetic field as we move toward the pole. Because the brightness 
is much lower at the equator than at the pole, a reasonable approximation is to consider 
the width constant with typically AR = 2 — 5 x 10 -3 (bottom-left panel) or even a slightly 
increasing width from the pole to the equator (bottom-right panel). A slightly increasing 
width will not change the overall morphology. It still produces a thin rim whose projected 
width is decreasing starting from the pole to the equator. In Figure [TT], we have assumed 
that the emissivity contrast at the BW between the pole and the equator is always a factor 
100. In fact, because the contrast between the pole and equator is larger in the X-rays 
than in the radio, we will obtain a similar morphology but with the size of the caps slightly 
reduced in X-ray. In other words, the X-ray synchrotron emitting rims become geometrically 
thinner and of smallest azimuthal extent as we go to higher energy. All these effects would 
be however slightly attenuated if we had included the instrumental effects (particularly the 
PSF), but overall this is consistent with what is observed in SN 1006. 



6.2. CR acceleration in a partially ionized medium 



When trying to delineate the BW, we found that the X-ray synchrotron-emitting rims 
and the diffuse filament of Ha emission were coincident over some little azimuthal range, 
notably for instance at the edge of the eastern cap where the synchrotron emission starts to 
turn on (see Fig. [T|). SN 1006 is one of the rare remnants where this characteristic is clearly 
observed, although this likely happens in the Tycho SNR too. 

In general, Ha emission from non-radiative shocks is believed to arise when the blast 
wave encounters partially ionized gas. This points to the presence of neutral atoms in the 
ambient medium around SN 1006 (jGhavamian et al.ll2002l ). Theoretical studies have shown 
that in such medium, the scattering Alfv en waves should be dam ped, henceforth quenching 
the acceleration of high energy particles (jO'C Drury et al.lll996l ). In fact, recent theoretical 
studies suggest that the m agnetic field could be turbulently amplified even in such medium 
f lBykov fc Toptyginll2005l ). The fact that we observe the Ha and synchrotron X-rays to be 
coincident gives support to this statement. Other observations of Ha emission coexisting 
with synchrotron X-rays would be very interesting in that regard. 
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6.3. Comparison between SN 1006 and Tycho 

Like SN 1006, the Tycho SNR is believed to be the remnant of a thermonuclear explosion 
and, as such, is expected to develop an approximately exponential ejecta density profile and 
to evolve in a uniform ambient medium. Measurements of the gap between BW an d CD 



have been done in the Tycho SNR using Chandra observations fjWarren et al.ll2005l ). We 
comment on the Tycho results in the light of our results. 

Because the thin and bright X-ray synchrotron-emitting rims concentrated at the BW 
are observed all around Tycho, in contrast to SN 1006, the BWs location can be easily 
determined. In addition, because the X-ray emission from the shocked ejecta is dominant 
in Tycho, in contrast to SN 1006, it is also fairly easy to determine the CD's location. The 
observed azimuthally-ayerage d ratio of radii between the BW and CD was found to be ~ 1.04 



in Tycho (IWarren et al.ll2005l ). After correcting, approximately, for effects due to projection 
of the highly structured CD, the ratio of radii became ~ 1.075. This is still a lot smaller than 
the (unprojected) ratio of ~ 1.25 derived using standard one-dimensional hydrodynamical 
simulations with no CR acceleration, at the current age of Tycho (430 yrs). It was then 
concluded that efficient CR ion acceleration was occurring around nearly the entire BW of 
Tycho. 



More recently, ICassam-Chenai et al.l (120071 ) applied a CR-modified self-similar hydrody- 



namic model to the radio and X-ray observations of the Tycho SNR using the observed ratio 
of BW to CD radii as an important input constraint. Using reasonable values for the hydro- 
dynamical parameters (ambient medium density, SN explosion energy, ejected mass), DSA 
parameters (injection efficiency, magnetic field, diffusion coefficient, electron-to-proton ratio 
at relativistic energies), and the distance, a good description of the observational properties at 
the blast wave of Tycho was obtained. This detailed study provides a self-consistent/coherent 
picture of efficient CR ion acceleration. 

Regarding SN 1006, we have some good news and some bad (compared to the case of 
Tycho, where CR-modified models seem to describe the data well). The good news is that 
there is a significant variation in the ratio of BW to CD radii as a function of azimuthal angle 
that is correlated with the varying intensity of the X-ray and radio synchrotron emission. 
Where the gap is large the synchrotron emission is faint (or not detected) and where the gap 
is small the emission is bright. This is consistent with an azimuthal variation in the effects 
of CR-modification to the remnant dynamics. The bad news is that in regions of presumably 
efficient acceleration (i.e., the bright rims), the ratio of BW to CD radii approaches small 
values, very close to 1. The puzzle is that our CR-modified hydrodynamical models are 
unable to produce ratio values less than 1.06, while the specific scenarios presented above 
produce ratio values closer to 1.12. In order to make progress in our study we decided to 
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introduce an ad-hoc rescaling of the ratio of radii by decreasing the modeled values at all 
angles by a constant factor of ~ 12%. This procedure works remarkably well and suggests 
to us that something else, unrelated to the CR acceleration process (such as 3-D projection 
effects and ejecta clumping), may be responsible for the smaller than expected gap between 
the BW and CD in SN 1006. 



7. Conclusion 

Using a combination of radio, optical (Ha) and X-ray images, we have located the 
positions of the BW and CD radii along the southeastern sector of SN 1006. We found that 
with increasing azimuthal angles the ratio of radii between the BW and CD increases from 
values near unity (within the northeastern synchrotron-emitting cap) to a maximum of about 
1.10 before falling again to values near unity (in the southwestern synchrotron-emitting cap). 
These variations reflect changes in the compressibility of the plasma attributed to variations 
in the efficiency of the BW at accelerating CR ions and give strong support to the overall 
picture that SNR shocks produce some fraction of Galactic CRs. However, at the present time 
we do not have a detailed explanation for the apparent overall smallness of the measured 
ratios of radii in SN 1006: the minimum value predicted by our CR-modified self-similar 
dynamical models is 1.06. In this study we simply rescaled the ratio of radii by a constant 
factor (independent of azimuthal angle) of ~ 12% with the expectation that some process, 
other than CR acceleration itself, was responsible for driving the edge of the ejecta closer to 
the BW all along the rim of SN 1006. The lack of a definitive astrophysical explanation for 
this discrepancy is a significant factor limiting our ability to understand the CR acceleration 
process in SN 1006. Further research into the effects of hydrodynamical instabilities at 
the CD and ejecta clumping (two possible explanations for the small BW/CD ratio) would 
incidentally provide new insights into the acceleration process. 

In addition to the azimuthal variations of the ratio of radii between the BW and CD 
we also interpreted the variations of the synchrotron flux (at various frequencies) at the BW 
using CR-modified hydrodynamic models. We assumed different azimuthal profiles for the 
injection rate of particles in the acceleration process, magnetic field and turbulence level. 
We found the observations to be consistent with a model in which these quantities are all 
azimuthally varying, being the largest in the brightest regions. Overall this is consistent 
with the picture of diffusive shock acceleration. In terms of morphology, we found that our 
model was generally consistent with the observed properties of SN 1006, i.e., a bright and 
geometrically-thin synchrotron-emitting rim at the poles and very faint synchrotron emission 
at the equator and in the interior. In the model, the X-ray synchrotron-emitting rims are 
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geometrically thinner and of smallest azimuthal extent than the radio rims, which is in broad 
agreement with observations. This is because the most energetic electrons accelerated at the 
blast wave lose energy efficiently in the amplified post-shock magnetic field. Based on this 
picture, it would be worth trying to measure the azimuthal variation of the magnetic field 
strength (from the X-ray rim widths) and thereby gain further insight into the amplification 
process. 
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Let r be the distance to the center O of a sphere of radius R s , (J) the azimuthal angle 
in the (Ox,Oz) plane with < <p < 2ir, and 9 the latitude with —tt/2 < 9 < it/ 2. The 
projection of a spherical emissivity distribution £(r, 9, 0) onto the plane (Ox, Oy) results in 
a two-dimensional brightness distribution: 



where I 2 = x 2 s — x 2 , and < x < x s = R s cos 9 and < y < R s . At any point of coordinate 
(x, y, z), the emissivity is £(x, y, z) = £(r, 9, 4>) where r 2 = x 2 + y 2 + z 2 and 9 = arcsin (y/r) 
and (f> = axccos(x / \/ x 2 + z 2 ). When the emissivity has a symmetry of revolution around the 
Oy axis (e.g., polar caps), there is no dependency on 0. 

Let us consider a sphere where the emissivity £ is radially decreasing from the maximum 
£max at the surface of the sphere with a characteristic width AR. In the case of an exponential 
decrease and cylindrical symmetry, we have: 



A. Projection along the line-of-sight in the case of polar caps 




(Al) 




(A2) 
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In the case of polar caps, £ max will be larger at the poles (8 = n/2 and 9 = — vr/2). 

In Figure [TTJ we show the projected morphology obtained with the emissivity spatial 
distribution given by Eq. (|A2[) in the range < 6 < n/2 assuming different values and 
angular dependencies for the width AR = AR/R S . In those plots, a squared-sine variation 
was assumed for log 10 (£ max ), having the value 0.0 at the pole and —2.0 at the equator. 
Decreasing the value at the equator produces similar rims but with a lower azimuthal extent 
(i.e., smaller polar caps). 
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Fig. 1. — Three-color Chandra image of SN 1006: 0.5-0.8 keV (red), 0.8-2 keV (green) and 2- 
4.5 keV (blue). The regions dominated by synchrotron emission from high-energy relativistic 
electrons appears in white and are naturally associated with the BW. The regions dominated 
by the thermal emission from the shocked gas (mostly oxygen) appears in red and are most 
likely associated with the ejecta. Note that the oxygen is found slightly ahead of the BW in 
some places (e.g., east and south). Point sources have been removed here so that only the 
diffuse emission from the remnant is visible. Images were background subtracted, corrected 
for vignetting and slightly smoothed. The intensity scaling is square-root with a maximum 
fixed at 5.0 x 10~ 4 , 5.0 x 10~ 4 and 2.0 x 10~ 4 ph/cm 2 /s/arcmin 2 for the red, green and blue 
bands, respectively. 
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Fig. 2. — Radio image of SN 1006 at 1.5 GHz. The intensity scaling is square-root with 
a maximum fixed at 1.0 mJy/beam. The bright elongated source in the western rim is a 
background radio galaxy. 




Fig. 3.- 



Ha image of SN 1006. 



Fig. 4. — Images of SN 1006 in several energy bands. Top-left: Radio image at 1.5 GHz. The 
maximum intensity was fixed at 0.6 mJy/beam. Top-right: Ha image with the contour errors 
associated with the shock front seen in the southeast (red lines). Bottom-left: Chandra X- 
ray image in the oxygen band (0.5-0.8 keV). The maximum intensity was fixed at 5.0 x 10~ 4 
ph/cm 2 /s/arcmin 2 . Bottom-right: Chandra image in the mid-energy X-ray band (0.8-2.0 
keV). The maximum intensity was fixed at 6.0 x 10~ 4 ph/cm 2 /s/arcmin 2 . Both X-ray images 
were background subtracted, corrected for vignetting and slightly smoothed. In the radio and 
X-ray images, we show the contour derived from the Ha emission (red lines). This contour 
lies within the contour errors shown in the Ha image. In all the images, the intensity scaling 
is linear. 
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Fig. 5. — Azimuthal variation of the outer Ha emission (black lines) and X-ray emission 
as obtained in several X-ray energy bands: 0.5-0.8 keV (red lines), 0.8-2 keV (green lines) 
and 2-4.5 keV (blue lines). Regions where the X-ray synchrotron emission is dominant are 
indicated by cross-hatched yellow lines. In those regions, the fingers (stars) indicate the 
presence of shocked ejecta found at or even slightly ahead of the BW. The X-ray contours 
correspond to places where the brightness becomes larger than 1.5 x 10 -5 ph/cm 2 /s/arcmin 2 
in the 0.5 — 0.8 keV band, 0.6 x 10~ 5 ph/cm 2 /s/arcmin 2 in the 0.8 — 2 keV band and 
0.4 x 10~ 5 ph/cm 2 /s/arcmin 2 in the 2 — 4.5 keV band. Note that the radii measurements 
are quite sensitive to the contour values, but the range of variations is contained between 
the red and blue lines. 
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Fig. 6. — Ratio of the outer Ha radius (black lines in Fig. [5]) to the 0.5-0.8 keV X-ray 
radius (red lines in Fig. [5]). Regions where the X-ray synchrotron emission is dominant 
are indicated by cross-hatched yellow lines. Error bars include uncertainties on the radius 
derived from the (expanded) Ha image and that done in the low-energy X-rays choosing 
different contour values (2.5 ± 1.0 x 10~ 5 ph/cm 2 /s/arcmin 2 ). The large error bars between 
200° — 205° are due to a low exposure there. 
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Fig. 7. — Top: Azimuthal variations of the brightness at the BW of SN 1006 at several 
frequencies/energies: 1.5 GHz (black lines), 0.5-0.8 keV (red lines), 0.8-2 keV (green lines) 
and 2-4.5 keV (blue lines). The brightness in each azimuthal bin was calculated between the 
BW radius and the radius 30" behind. For comparison, we show the level of the radio rms 
noise level and X-ray background located 1' outside the BW and extracted from a 30"-wide 
box (dashed lines). Bottom: Same but with the curves rescaled roughly to the same level in 
the brightest regions and shown with a log scale. 
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Fig. 8. — Ratio of radii between the BW and CD, R^w/Rcb, & s obtained from 1-D numerical 
hydrodynamical simulations (assuming no CR acceleration at the BW) and calculated for a 
large range of ambient density, no, for a fixed SNR age (dashed lines: 400 yrs; solid lines: 
1000 yrs). We show curves obtained with a 1-D exponential (EXP) ejecta profile with a 
kinetic energy of the explosion E51 = -E/10 51 ergs = 1.0 (red lines), E$i = 0.8 (blue lines), 
E 5 i = 1.4 (green lines) and a 1-D powerlaw (PL7) ejecta profile (n = 7) with E 5 i = 1.0 
(yellow lines). The cross-hatched domains correspond to the intersection of the range of ratio 
of radii observed in SN 1006 and Tycho (< 1.10, vertical solid lines) and range of ambient 
medium density typically encountered in their vicinity (horizontal dotted lines). The 1-D 
hydrodynamical simulations are clearly unable to predict a ratio of radii as small as the 
one observed in both SNRs. For completeness, we indicate by a filled circle the value of 
the ambient density at which the simulated SNR radius matches the one observed for each 
model. An angular BW radius of 14.6' (resp. 256") and a distance to the SNR of 2.0 kpc 
(2.8 kpc) were assumed for SN 1006 (resp. Tycho), resulting in a physical BW radius of 
~ 8.5 pc (resp. ~ 3.5 pc). Larger values for the physical radius (due for instance to a larger 
distance to the remnant) would correspond to lower ambient densities. 
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Fig. 9. — Various models assuming different azimuthal profiles for the injection rate, rj m j 
{top), and far upstream magnetic field, B n (second, solid lines). Each model predicts az- 
imuthal variations of the ratio of radii between the BW and CD, Rbw/Rcd (third), and 
synchrotron brightness (projected along the line-of-sight) at different frequencies (bottom). 
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Fig. 10. — Left column: azimuthal profiles assumed for the injection efficiency, ?7i n j (panel 
a), far upstream magnetic field, B u (panel b) and diffusion coefficient normalized to the 
Bohm-limit value, ko (panel c). Middle column: model profiles of the overall density (r to t) 
and magnetic field (r B ) compression ratios (panel d), immediate post-shock value of the 
magnetic field, B d (panel e) and maximum energies (panel /). Right column: predicted 
(thick lines) and observed (thin histograms) azimuthal profiles for the ratio of radii between 
the BW and CD, Rbw/Rcd (panel g, solid lines) and line-of-sight projected synchrotron 
brightness at several frequencies (panels h and i). 



Fig. 11. — Projected morphology predicted in the case of polar caps (see Appendix [AT) 
with AR = 5 x 1(T 2 (top-left), AR = 10" 2 (top-right), AR = 5 x 10" 3 (bottom-left), 
AR linearly varying from 2 x 10 -3 at the pole to 8 x 10 -3 at the equator (bottom-right). 
The radio and the high-energy X-ray morphology predicted by the model shown in Figure 
[TU1 correspond roughly here to the top-right and bottom-right panels, respectively. In all 
panels, the intensity scaling is square-root. 



